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O I A search for new bosons possessing couplings to lepton-quark pairs is performed in the 

HI experiment at HERA using 1994 to 1997 data corresponding to an integrated luminos- 
^ . ity of 37pb First generation leptoquarks (LQs) are searched in very high neutral 

^ \ (NC) and charged (CC) current data samples. The measurements are compared to Stan- 

><^ ■ dard Model (SM) expectations from deep-inelastic scattering (DIS). A deviation in the 

spectrum previously observed in the 1994 to 1996 dataset at Q^-^ 15000 GeV^ remains, 
OhI though with less significance. This deviation corresponded to a clustering in the invariant 

^ ' mass spectrum at M ~ 200 GeV which is not observed with the 1997 dataset alone. The 

NC DIS data is used to constrain the Yukawa couplings A of first generation scalar and 
vector LQs in the Buchmuller-Riickl-Wyler effective model. Scalar LQs are excluded for 
^ ' masses up to 275 GeV for a coupling of electromagnetic strength, A = 0.3. A sensitivity to 

■ coupling values < 1 is established for masses up to 400 GeV for any LQ type. The NC and 

CC DIS data are combined to constrain A for arbitrary branching ratios of the LQ into eq 
in a generic model. For a decay branching ratio into e'^u pairs as small as 10%, LQ masses 
up to 260 GeV are ruled out for A = 0.3. LQs possessing couplings to mixed fermion 
generations, which could lead to signals of lepton flavor violation (LFV), aie searched in 
events with a high transverse momentum /i or r. No ^ + X or t + X event candidate is 
found that is compatible with LQ kinematics. Constraints are set on the Yukawa coupUng 
involving the ^ and r lepton in a yet unexplored mass range. 
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1 Introduction 



The ep collider HERA offers the unique possibility to search for s-channel production of new 
particles which couple to lepton-parton pairs. Examples are leptoquark (LQ) colour triplet 
bosons which appear naturally in various unifying theories beyond the Standard Model (SM) 
such as Grand Unified Theories [[l|] and Superstring inspired Eq models [|^], and in some Com- 
positeness 0] and Technicolour [Q] models. Leptoquarks could be singly produced by the fusion 
of the 27.5 GeV initial state lepton with a quark of the 820 GeV incoming proton, with masses 
up to the kinematic limit of ^/s^ — 300 GeV. 

The interest in such new bosons has been considerably renewed recently following the ob- 
servation by the HI ^ and ZEUS [p experiments of a possible excess of events at very high 
masses and squared momentum transfer Q^, above expectations from SM neutral current (NC) 
and charged current (CC) deep-inelastic scattering (DIS). These early results were based on 
data samples collected from 1994 to 1996. Of particular interest was the apparent "clustering" 
of outstanding NC events at masses around 200 GeV observed in HI which has motivated con- 
siderable work on LQ kinematics constraints and phenomenology [|^], extending beyond 
the original effective model of Buchmiiller-Riickl-Wyler (BRW) [^. 

In this paper, LQs are searched using all available e^p data collected in HI from 1994 to 
1997. Inclusive single and double differential DIS cross-sections obtained from a similar dataset 
are presented in a separate paper Here, firstly, NC and CC measurements at high are 
compared with SM expectations at detector level. Mass and angular distributions of NC- and 
CC-like events are then used to set constraints on first generation LQs. The search is then further 
extended to LQs possessing couplings to leptons of different generations. Such lepton flavor 
violating (LEV) LQs would lead to final states involving a second or third generation lepton. 

The total integrated luminosity C amounts to 37 pb \ an increase in statistics of a factor 
~ 2.6 compared to previous HI analysis at very high []^] and a factor ~ 13 compared to 
previous LQ searches at HERA [jlT], [l^, [Dp . 



2 The HI Detector 

A complete description of the HI detector can be found elsewhere M]. Here we introduce 



only the components relevant for the present analysis in which the final state of the processes 
involves either a charged leptonQ with high transverse energy or a large amount of hadronic 
transverse energy flow. 

Positron energies and angles are measured in a liquid argon (LAr) sampling calorimeter [ [T3| ] 
covering the polar angulaiQ range A° <0 < 154° and all azimuthal angles. The LAr calorimeter 
is divided in eight "wheels" along the beam axis, themselves subdivided in up to 8 modules with 
minimum inactive material ("cracks") in between. The modules consist of a lead/argon elec- 
tromagnetic section followed by a stainless steel/argon hadronic section. Their fine read-out 

'The analysis does not distinguish expHcitly between + and — charges. 

^The z axis is taken to be in the direction of the incident proton, the forward direction, and the origin of 
coordinates is the nominal ep interaction point. 
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granularity is optimized to provide approximately uniform segmentation in laboratory pseudo- 
rapidity and azimuthal angle 0. Electromagnetic shower energies are measured with a resolu- 
tion of a{E)/E ~ 12%/ ^E/ GeV © 1% and pion induced hadronic energies with a{E)/E ~ 
50%/ ^yE/ GeV © 2% after software energy weighting. These energy resolutions were mea- 
sured in test beams with electron energies up to 166 GeV [ [T6| , [T^] and pion energies up to 
205 GeV [|T7|]. The energy calibration was determined initially from test beam data with an un- 
certainty of 3% and 4% for electromagnetic and hadronic energies respectively. A new absolute 
energy scale calibration for positrons detected in the actual HI experiment has been recently 



established [ [T0| , [181 ] in situ by using the over-constrained kinematics of NC DIS, QED Comp- 
ton and e+e^ pair production from two-photons processes. A precision of 0.7% is reached in 
the LAr central barrel region 80° <9e< 145°, 1.5% in 40° < 6^ < 80° and 3.0% in the forward 
region 5° < < 40°. The precision on the hadronic energy scale was determined by requiring 
the balance of the transverse momenta of the positron and hadronic system in NC DIS events. 
This was performed using a method [ [T0| ] correcting the energy flow associated to jets by LAr 
wheel calibration constants. The hadronic energy scale is found to be understood at the 2% 
level when comparing to Monte Carlo expectation. This represents an improved understand- 
ing of both the electromagnetic and hadronic energy scales compared to [Q], made possible 
by the increase of statistics accumulated in 1997. All analyses described in the following rely 
on this updated calibration. The resolution on the polar angle of the positron measured from 
the electromagnetic shower in the calorimeter varies from ~ 2 mrad below 30° to < 5 mrad 
at larger angles. A lead/scintillating-fibre backward calorimeter [ |T9| ] extends the coverage[| at 
larger angles (153° <9< 178°). 

Located inside the calorimeters is the tracking system which is used here to determine the 
interaction vertex and provide charged track information relevant for lepton identification (see 
section The main components of this system are central drift and proportional chambers 
(25° <9< 155°), a forward track detector (7° < 6* < 25°) and a backward drift chamber^ The 
tracking chambers and calorimeters are surrounded by a superconducting solenoid providing a 
uniform field of 1.15 T parallel to the z axis within the detector volume. The instrumented iron 
return yoke surrounding this solenoid is used to measure leakage of hadronic showers and to 
recognize muons. The luminosity is determined from the rate of the Bethe-Heitler ep ep'j 
bremsstrahlung measured in a luminosity monitor. This consists of a positron tagger and a 
photon tagger located along the beam line, —33 m and —103 m respectively from the interaction 
point. 

For the acquisition of events we rely on the timing information from a time-of-flight system 
and on the LAr trigger system which provides a measurement of the energy flow using coarse 
trigger towers [jT5|]. 



3 Leptoquark Phenomenology and Models 

Leptoquark production at HERA can lead to final states similar to those of DIS physics at very 
high Q^. The basic DIS processes are illustrated in Fig. |l]a. Leptoquarks can be resonantly 

''The detectors in the backward region were upgraded in 1995 by the replacement of the lead/scintillator 
calorimeter ^Bj and a proportional chamber. 
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produced in the s-channel and exchanged in the w-channel as illustrated by the diagrams in 
Fig. [l]b,c. Here and in the following, whenever specified, the indices i and j of the couplings \ij 
at the LQ—lepton— quark vertices refer to the lepton i*^and quark j*^ generation respectively. 
Otherwise, A designates a coupling of LQs to first generation fermions. LQs coupling only 



Figure 1: Diagrams of (a) deep-inelastic scattering; (b) s-channel resonant production and (c) u- 
channel exchange of a leptoquark with fermion number F = 0. Diagrams involving a | F |= 2 lepto- 
quark are obtained from (b) and (c) by exchanging q and q. 

to first generation fermions (henceforth called first generation LQs) give e + q or u + q' final 
states leading to individual events indistinguishable from SM NC and CC DIS respectively. 
LQs with LFV couplings to second or third generation leptons (henceforth called LFV LQs) 
can participate in e~^p ^ /i"^ + jet + X or e~^p — > r+ + jet + X processes. Such exotic 
signatures are expected to be essentially background free for high transverse momentum of the 
observable lepton. 

In the s-channel, a LQ is produced at a mass M = ^/s^ where x is the momentum 
fraction of the proton carried by the struck quark. Over a large fraction of the mass range 
accessible at HERA and for a reasonable coupling strength, e.g. satisfying A^/47r < 1, the 
intrinsic decay width of a scalar (S) or vector (V) LQ of nominal mass M^q into a lepton and 
a quark is expected to be small. This width is calculated as Ts = {3/2)Tv = X^jM^q /16n 
which corresponds for example to r5 ^ 40 MeV for a scalar at Mlq = 200 GeV and Ajj = 
O.L In the narrow-width approximation (NWA), the resonant production cross-section ctnwa is 
proportional to A^g(x) where q(x) is the density of the struck parton in the incoming proton. 
However when approaching the kinematic limit where the values of q{x) are very small, the 
coupling strengths which can be probed with the actual integrated luminosities are too high 
for the NWA to be valid. The convolution of the steeply falling q(x) with the Breit-Wigner 
distribution of finite width characterizing the resonance leads to a strong distortion of the LQ 
mass peak, and the mass spectrum shows very large tails towards low values. As a result the LQ 
production cross-section cr^ in the s-channel for M^q approaching ^/s^ is considerably larger 
than (Tnwa- The deviation from (Tnwa is significant (typically > 10%) at Mlq > 250 GeV for a 
LQ produced via a valence quark (u, d), and already at Mlq > 200 GeV for a production via 




(a) 



(b) 



(c) 
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a sea quark (u, d) pTl]. The analysis presented in the following fully takes into account these 
effects originating from the finite LQ decay width. 

Scalar LQs produced in the s-channel decay isotropically in their rest frame leading to a flat 
da j dy spectrum where y = Q'^/s^pX = | (1 + cos^*) is the Bjorken scattering variable in 
DIS and 9* is the decay polar angle of the lepton relative to the incident proton in the LQ centre 
of mass frame. In contrast, events resulting from the production and decay of vector LQs would 
be distributed according to da / dy oc (1 — y)^. These y spectra (or in other words the specific 
angular distributions of the decay products) from scalar or vector LQ production are markedly 
different from the da / dy oc distribution expected at fixed x for the dominant t-channel 
photon exchange in neutral current DIS events Q. Hence, a LQ signal in the NC-like channel 
will be statistically most prominent at high y. 

The ?i-channel contribution scales with A^. It can compete with resonant production only 
for LQs with fermion number \ F \= 2 and at high couplings and LQ masses. For F = 
LQs, it is highly suppressed by less favorable parton densities as it proceeds via an exchange 
involving an antiquark from the proton. Scalar LQ exchange would lead to events distributed in 
y according to da / ~ {1 — yY while vector LQ exchange would lead to a flat y spectrum. 
However the events originating from w-channel LQ exchange would mainly be concentrated at 
mass values much lower than M^q. As such, the kinematic cuts used in this analysis to reduce 
the number of NC-like and CC-like events (see section ||) also drastically suppress a possible 
M-channel contribution. 

In approaching Mlq ~ ^/s^, the interference of the LQ s-channel production and u- 
channel exchange with SM boson exchange can no longer be neglected. This interference can 
be constructive or destructive^ depending on the LQ type. As will be seen in section ^ the set 
of cuts used in the present analysis focuses on a phase space region where the contribution of 
the interference is considerably reduced. 





Angular Spectrum 


Cross-Section 




y Shape 


A Dependence 




Scalar 


Vector 


M < 


M > 


s-channel 


flat 






X2 X2 


M-channel 




flat 




X2 X2 


Interference 






A2 



Table 1 : Main properties of the different LQ induced contributions at HERA to e+ + — > + qj. For 
LQs coupling to both eq and /„g pairs with /„ 7^ e, (3n denotes the branching ratio of the LQ into /„ + qj. 
The interference contribution only concerns processes with a first generation lepton in the final state. 

The experiments at HERA are also sensitive to LQs with Mlq ^ ^/s^. For first generation 
LQs, the interference between LQ induced and SM boson exchange processes (which scales 

'^At high momentum transfer, exchange is no longer negligible and contributes to less pronounced differ- 
ences in the y spectra between LQ signal and DIS background. 

^ The signs of the interference terms between SM gauge boson and F = LQ contributions given in the 
original BRW paper [||] were found to be incorrect. The correct signs as provided in the erratum to Ref. have 
been used here. 
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with A^) generally dominates in this mass range over the s-and w-channel contributions (both 
scaling with A^). For F = LQs, a nevertheless sizeable s-channel contribution originates from 
the convolution of the parton density with the low mass tail of the LQ Breit-Wigner resonance 
of finite width, provided that the coupling A is not too small. For example, for A = 1 and for the 
LQ labelled 5'i/2,l in the BRW model (see below), the s-channel contribution competes with 



the interference for M^q up to ~ 400 GeV, within the kinematic cuts used in section q4 



In the cases of M^q ^ y/s^, the propagator entering the LQ amplitudes can be con- 
tracted to a four-fermion interaction. One is left with a contact interaction mostly affecting the 
measured inclusive DIS spectrum through interference effects. Constraints on such four- 
fermion couplings translated into limits on M^q/X for first generation LQs will be discussed 
in a separate paper. For LFV LQs above y^s^, both the s- and w-channel contributions may 
in principle be important for large Yukawa coupling values. There, the LFV LQ cross-sections 
a{eqi InQj) (s-channel) and a{eqj InQi) (w-channel) only depend on Aii, A„j and Mlq via 



Some essential characteristics of the different LQ induced processes contributing at HERA 
are summarized in table [T]. 

The LQ searches will be discussed here either in the strict context of the BRW phenomeno- 
logical ansatz where the decay branching ratios are fixed by the model, or in the context of 
generic models allowing for arbitrary branching ratios. The BRW model considers all possible 
scalar (Sj) and vector (Vj) LQs of weak isospin / with dimensionless couplings X-J to lepton- 
quark pairs, where L or 7? is the chirality of the lepton. The general effective Lagrangian which 
is introduced obeys the symmetries of the SM. There are 10 different LQ isospin multiplets, 
with couplings to left or right handed fermions, among which there are 5 isospin families of 
scalar LQs. These are listed in table ^ We restrict the search to pure chiral couplings of the 
LQs given that deviations from lepton universality in helicity suppressed pseudoscalar meson 
decays have not been observed [^, This restriction to couplings with either left- (A^) or 
right-handed (A^) leptons (i.e. A^ ■ A^ ~ 0), affects only two scalar (Sq and 51/2) and two vec- 
tor (Vi/2 and Vo) LQs. We make use of the so-called Aachen nomenclature and classification 
scheme [|2^] and do not use specific symbols to label the anti-leptoquarks which are actually 
produced in e+p collisions. We make the simplifying assumptions that one of the LQ multiplets 
is produced dominantly and that the mass eigenstates within the LQ isospin doublets and triplets 
are degenerate in mass. 

For the determination of LQ signal detection efficiencies, we make use of the LEGO event 
generator p5| ] and of a complete simulation of the HI detector response. LEGO incorporates 
s- and M-channel LQ exchange processes depicted in Fig. [l|b,c. It takes into account initial 
state QED radiation in the coUinear approximation. The parton showers approach y26| J relying 



on the DGLAP p7| ] evolution equations is used to simulate QCD corrections in the initial and 
final states, and the kinematics at the decay vertex is properly corrected for effects of the parton 
shower masses. The non-perturbative part of the hadronization is simulated using string frag- 



mentation [pq]. The Mandelstam variable s characterizing the eg Iq subprocess defines the 
scale at which the parton density is evaluated as well as the maximum virtuality of parton show- 
ers. The LQ signal cross-sections are calculated using the full matrix elements given in Ref. [|^] 
and taking into account the contributions from the s- and w-channels as well as the interference 



8 



F = -2 


Prod./Decay 


(3e 


F = 


Prod./Decay 




Scalar Leptoquarks 


'/'So 




1/2 


'/'S,/2 


(i^UR e+M 


1 




e~I^UL e+n 


1 




e[UL e+?i 


1 


'/'So 




1 


'/'Si/2 


e^^di — > e+(i 


1 




e^dji e^d 


1 


'/'~Sl/2 


e^ciij e^d 


1 




e%UR e+M 


1/2 








Vector Leptoquarks 


'/'Vy2 


e^rf/j e^d 


1 


'/'Vo 


e~ldfi e^d 


1 




e^Ji — > e+rf 


1 




e^di — i> e^d 


1/2 


'/'Vy2 




1 


'/'% 


etuR e+M 


1 


'/'Vl/2 




1 




e~^UL e~^u 


1 










e~^dL — >■ e~^d 


1/2 



Table 2: Leptoquark isospin families in the Buchmiiller-Riickl-Wyler model. For each leptoquark, the 
superscript corresponds to its electric charge, while the subscript denotes its weak isospin. For simplicity, 
the leptoquarks are conventionally indexed with the chirality of the incoming electron which could allow 
their production in e~p collisions, e.g. the So will be denoted by So^r (see text). /3e denotes the branching 
ratio of the LQ into + q. 



with SM boson exchange. The resulting cross-sections are further corrected to account for next- 
to-leading order (NLO) QCD effects making use of multiplicative /^-factors pB| ] in a procedure 
described in detail in section These NLO QCD corrections which depend on the LQ signal 
shape expected for a given Mlq and A are typically of 0(10%). For the parton densities, use 
is made of the recent Martin-Roberts-Stirling-Thome MRST [ p9[ ] parametrization which better 
describes existing measurements constraining the sea quark densities in the proton PC|, |3T[]. 

The theoretical uncertainty on the signal cross-section originating mainly from contributions 
of parton density distributions extracted from "QCD fits" and the value of the strong coupling 
constant «5 is treated as a systematic error. This uncertainty is ~ 7% for leptoquarks coupling to 
e'^u, and varies between ~ 7% at low LQ masses up to ~ 30% around 250 GeV for leptoquarks 
coupling to e^d. Above 250 GeV, for coupling values corresponding to the expected sensitivity, 
the small but finite width of the resonance results in the fact that mainly relatively low x partons 
are involved in the LQ production. Hence the uncertainty on the signal cross-section decreases 
for F = LQs to ~ 7% between 250 GeV and the kinematic limit. For \ F \= 2 leptoquarks, 
this uncertainty ranges from ~ 10% at low masses and reaches ~ 40% around 200 GeV, and 
then goes down to ~ 15% at the kinematic limit. Moreover, choosing alternatively or the 
square of the transverse momentum of the final state lepton instead of s as the hard scale at 
which the parton distributions are estimated yields an additional uncertainty of ±7% on the 
signal cross-section. 
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4 Deep Inelastic Scattering and Other Background Sources 



The calculation of the SM expectation for NC and CC DIS ep scattering is performed using the 
parton model in the approximation of single •y/Z and W boson exchange, and relies on a de- 
scription of the proton in terms of scale dependent structure functions. The structure functions 
are expressed in terms of parton densities and are taken here from the MRST parametrization 
which includes constraints from HERA data up to Q'^ = 5000 GeV^ [521 530 . The parton densi- 



ties are evolved to the high Q"^ domain relevant for this analysis using the next-to-leading order 
DGLAP equations. The Monte Carlo event generator DJANGO [B^ which follows such an 



approach is used for the comparison with data. This generator includes the QED first order ra- 
diative corrections and a modelling of the emission of QCD radiation via ARIADNE pB]. 



The ARIADNE generator makes use of the Colour Dipole Model [ p7| | to simulate QCD radia- 
tion to all orders and string fragmentation to generate the hadronic final state. 

The contributions from all background processes which could give rise to events with true 
or misidentified isolated leptons at high transverse energy or to events with a large missing 
transverse momentum have been evaluated. In particular, direct and resolved photoproduction 
processes were modelled using the PYTHIA generator [^Sp. It is based on leading order QCD 
matrix elements and includes initial and final state parton showers calculated in the leading log- 
arithm approximation, and string fragmentation. The renormalization and factorization scales 
were both set to P^, Pt being the transverse momentum of the jets emerging out of the hard 
subprocess. The GRV (G) leading order parton densities in the proton (photon) have been 
used pP|]. The production of electroweak vector bosons and was modelled using the 



EPVEC event generator. Contributions from two-photon processes where one 7 originates 
from the proton were also considered and estimated using the LPAIR [ ^T] ] event generator. A 
complete Monte Carlo simulation of the HI detector response has been performed for all back- 
ground processes. 

The following experimental errors are propagated as systematic errors on the mean SM 
expectations : 

• the uncertainty on the integrated luminosity (±1.5%); 

• the uncertainty on the absolute calibration of the calorimeters for electromagnetic ener- 
gies, ranging between ±0.7% in the central LAr wheels to ±3% in the forward region of 
the LAr calorimeter (see section 

• the uncertainty on the calibration of the calorimeters for hadronic showers of ±2% (see 
section 

In addition, a 7% theoretical uncertainty on the predicted NC DIS cross-section originates 
mainly from the lack of knowledge on the proton structure (see detailed discussion in |Q]) and, 
to a lesser extent, from the higher order QED corrections. For CC DIS processes which are 
mainly induced by d quarks, this uncertainty varies with Q"^ and ranges between 7% and ~ 20% 
at the highest considered here. All analyses described in the following sections have been 
repeated with an independent shift of the central values by ±1 standard deviation of each of 
the experimental and theoretical sources of errors. The overall systematic error of the SM 
prediction is determined as the quadratic sum of the resulting errors and of the statistical error 
on the Monte Carlo simulation. 
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5 Event Selection and Comparison with Standard Model Ex- 
pectation 



The search for first generation LQs relies essentially on an inclusive NC selection requiring an 
identified positron at high transverse energy, and an inclusive CC selection requiring a large 
missing transverse momentum. The selection of NC- and CC-like events will be described 



in subsections |5n| and respectively where the resulting samples will be compared to SM 



expectations. For LQs possessing couplings to mixed fermion generations, the selection of 



e'^p —* fi^ + q+X ande~^p r+ + g+X candidates which will be discussed in subsections p3 
and ^]4] respectively, requires essentially an identified /i or r lepton together with a large amount 
of hadronic transverse energy. 

In common for all channels analysed, the events must have been accepted by a LAr trigger 
asking either for an electromagnetic cluster, for a large transverse energy in the central part of 
the calorimeter, or for a large imbalance in the transverse energy flow. The rejection of back- 
ground from cosmic rays and from "halo" muons associated with the proton beam mainly relies 
on constraints on the event timing relative to the nominal time of the beam bunch crossings. 
Beam-wall and beam-residual gas interactions are furthermore suppressed by requiring a pri- 
mary vertex in the range \ z — z \ < AO cm where z varies within ±5 cm around z = depending 
on the HERA beam settings. 

In what follows, unless explicitly stated otherwise, the energy flow summations run over all 
energy deposits i in the calorimeters (apart from the electron and photon taggers). Thus, the 

missing transverse momentum PT,miss is obtained as PT,miss = \J Ex,if + (XI ^y,if with 
Ex^i = Ei sin 6i cos (pi and Ey^i = Ei sin 6i sin 0j. The momentum balance with respect to the 
incident positron is obtained as ^ (i? — Pz) = Yl i^i ~ Ez,i) with E^^i = Ei cos 9i. 



5.1 Neutral Current Deep-Inelastic-like Signatures 
5.1.1 Event selection and kinematics 

The selection of NC DIS-like events uses mainly calorimetric information for electron finding 
and energy-momentum conservation requirements, with selection cuts similar to those consid- 
ered in previous high Q'^ analysis @]: 

1. an isolated positron with Ex^e > 15 GeV (-Er.e = -EeSin6'e), found within the po- 
lar angular range 5° < < 145°; the positron energy cluster should contain more 
than 98% of the LAr energy found within a pseudorapidity-azimuthal cone of opening 
a/ (Ar^e)^ + (A0e)^ = 0.25 where r^e = — In tan ^; at least one charged track is required 
within the positron isolation cone; 

2. a total transverse momentum balance Pr.miss/ — 4vGeV ; 

3. a limited reconstructed momentum loss in the direction of the incident positron such that 
40 GeV <Y.iE - P,) < 70 GeV. 
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The identification of positron induced showers relies on the detailed knowledge of the expected 



lateral and longitudinal shower properties [ [161 , [ISl l. The efficiency for the detection of positrons 



exceeds 90% everywhere within the acceptance cuts, the main losses being due to showers de- 
veloping through the inactive material between calorimeter modules. The cut (2) makes possible 
a very efficient NC DIS selection up to the highest by taking into account the natural depen- 
dence of the calorimetric energy resolution, Ex^e being used as an estimate of the scale relevant 
for the actual PT,miss measurement. The cut (3) retains more than 90% of NC DIS events and 
exploits the fact that by energy-momentum conservation, the ^{E — P^) distribution for NC 
DIS events is peaked at 2E^, where E^ is the positron beam energy. It rejects events where a 
very hard coUinear 7 is emitted by the initial state positron. To ensure a good control of the 
positron identification performances, a fiducial cut is applied requiring: 

4. an azimuthal impact of the track associated to the positron at\ (p^ — (j)crack I > 1° from the 
nearest projective </> crack in the transverse plane. 

The DIS Lorentz invariants Q'^, y and M are determined using only the measurement of the 
"scattered" positron energy and angle as soon as | 0e — (pcrack |> 2°, such that the measurement 
of -Eg is reliable : 

M. = ,/^, Ql- ■ E,-EeCose, 



This method will henceforth be called the electron method (e-method). In ~ 4% of the accep- 
tance corresponding to the range 1° <| 0e — (j^crack \< 2° where calorimetry measurements of 
positrons deteriorate, the reconstructed positron energy is corrected to the value given by the 
double angle method [H2n : 

2E° 1 



E. 



2a 



ae + dh sin 6e 



usmg 



a, = tan(ee/2) = ^ and ah = tan(0ft/2) - ^'^^^ 



\2 



where the summations run over all energy deposits of the hadronic final state. 

In the following analysis, the comparison with SM expectation is restricted to the kinematic 
range > 2500 GeV^ and 0.1 < y < 0.9. The resolution in Me degrades with decreasing ye 
{5 Me /Me oc l/^/e) and so the low y domain is excluded. Excluding the high y values avoids 
the region where migrations effects due to QED radiation in the initial state are largest for the 
e-method. In the kinematic range considered and given cuts (1) to (4), the NC trigger efficiency 
exceeds 98% and is consistent with 100% to within experimental error. 

The y < 0.9 restriction also suppresses the photoproduction background where e.g. a jet 
has been misidentified as an electron. Following [^, any possibly remaining non-DIS contam- 
ination coming for example from 77 or QED Compton processes as well as the background 
from misidentified low NC DIS are further suppressed through a minimal set of specific 



cuts [18]. Among these a prominent one against multi-lepton final states is the requirement of 
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at least one reconstructed jet with Ex jet > 7 GeV found using a cone algorithm in the labora- 
tory frame, with a radius a/ (Ar/)^ + (A^)^ = 1. The jet should be in the polar angular range 
7° < Ojet < 145°, and at least 5% of its energy should be deposited in the hadronic section of 
the LAr calorimeter. Against photoproduction and low Q"^ NC DIS, it is also required that there 
be less than 7.5 GeV in the backward calorimeter. The specific reduction of the background 
contamination induces only small (< 5%) efficiency loss for high NC DIS-like processes. 
The remaining contamination is estimated to be below 0.15% and is henceforth neglected. 

Leptoquark signal selection efficiencies are determined over a coupling-mass grid with steps 
in mass of 25 GeV and for coupling values corresponding roughly to the expected sensitivity 
to properly take into account the effect of the intrinsic finite width of the searched resonance. 
Detailed Monte Carlo simulation of about 500 events per point on the grid is performed followed 
by the application of the full analysis chain. The above set of cuts ensures a typical selection 
efficiency for LQ — > e + g events which varies between ~ 40% and ~ 75% for LQ masses 
ranging in 75 to 250 GeV. 

Applying all the above NC selection criteria, 1298 DIS event candidates are accepted which 
is in good agreement with the expectation of 1243 ± 95 events from standard NC DIS. 

5.1.2 Comparison with Standard Model expectation 

Figure ^ shows the distribution of the NC candidates in the ye - Me kinematic plane. In such a 
plane, the signal of a 200 GeV scalar LQ with F = could manifest as illustrated in Fig. for 
coupling values corresponding to the expected sensitivity. Compared to other commonly used 
kinematic methods for NC DIS at HERA []TD[], the e-method provides the best peak resolution 
(truncated Gaussian fit) in mass at high y, where a LQ signal would be most prominent. This 
resolution on Me varies within ~ 3 — 6 GeV for LQ masses ranging between 100 and 250 GeV. 
It should be noted, as can be inferred from Fig. |^, that the e-method underestimates on average 
the true LQ mass by 2% due to migrations caused by final state QCD radiation. 

A differential analysis in the ye- Me plane of the very high events from the 1994 to 1996 
datasets [j^] had revealed a noteworthy excess of NC DIS-like events at Ql > 15000 GeV^ or for 
ye > 0.4 at masses Me — 200 GeV. A comparison of the measured M, y and Q"^ spectra with 
standard DIS model expectations can now be re-examined with higher statistics. 

We first consider the Me and ye information. The projected Mg and ye spectra are shown 
in Fig. |3| in several kinematic domains. Figures ^ and ^ show the projected Mg and ye dis- 
tributions of the NC DIS-like selected events at "moderate" (2500 < Ql < 15000 GeV^) 
and Fig. |c and |d at "very high" {Ql > 15000 GeV^). The distributions of the measured 
data are well reproduced by standard DIS predictions in the low range. At high the data 
slightly exceed the NC DIS expectation at Mg ~ 200 GeV as can be seen in Fig. |3|c. Moreover, 
Fig. 1^ shows that the excess of observed events is more prominent at high ye, so that at high 
Me and large ye the experiment tends to exceed the SM expectation. 

Figure ^ shows the measured and expected Mg distributions for a minimum ye value of 
ymin = 0.4. An excess of events over the NC DIS expectation at high mass (~ 200 GeV) is 
still visible. In the mass range 200 GeV ± AM /2 with AM = 25 GeV, Nobs = 8 events are 
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Figure 2: Kinematics in the ye - plane of (a) the selected NC DIS candidate events from HI data 
(two isocurves at = 2500 and 15000 GeV^ are plotted as full lines); (b) a scalar F = leptoquark of 
mass Mlq = 200 GeV decaying into e + q, for a coupling A = 0.05. 



observed for an expectation of Ndis = 2.87 ± 0.48. The mean mass value of these 8 events as 
determined with the e-method (Me) = 202.5 ± 7.0 GeV (RMS) agrees within 1.4% with the 
one obtained from the invariant mass of the final e-jet pairs. Of the observed events, 5 originate 
from the 1994 to 1996 data (40.3% of C) and 3 from the 1997 data (59.7% of £). It should be 
emphasized here that Nohs and Ndis are quoted for the same AA/ — Ay region where the most 
significant excess was observed in the original analysis of the 1994 to 1996 data [|5|] despite 
the fact that the individual events are slightly (within originally estimated systematic errors) 
displaced in the M — y plane. In this domain, 7 events were reported in [^. These events are 
measured here at values on average 2.4% higher due to the new in situ calibration of the 
electromagnetic section of the LAr calorimeter, and thus, one event has now migrated outside 
this AM domain. The estimated mass of one of the other 6 events in which the positron lies 
within less than 2° from the closest (p crack was and remains measured outside this mass region 
when using the double angle method in contrast to the e-method used in [^. It was explicitly 
checked that repeating the 1994 to 1996 analysis procedures of ^ but using this new calibration 
leaves the statistical significance and the physics messages of [|]] unchanged. 

At large mass Me > 180 GeV and for ye > 0.4, we observe in the 1994 to 1996 data 
Nobs = 7 in slight excess of the expectation of N^js = 2.21 ± 0.33 while in the 1997 data alone 
^obs = 4 events are observed, in good agreement with the expectation of N^js = 3.27 ± 0.49. 



14 



c 


400 


!5 




in 


300 


c 




— ' 


200 








100 




n 
u 


c 




!5 


10 


in 




c 


7.5 






> 


5 




2.5 








2500 < Q- < 15000 GeV^ 



2500 < Qc < 15000 GeV^ 



!5 



c 

0) 




100 150 200 250 300 
(GeV) 

> 15000 GeV^ 




100 150 200 250 300 
(GeV) 

y >0.4 



10 



a; 10 




100 150 200 250 
(GeV) 




c 

in 
c 
> 

0) 



0.2 0.4 0.6 0.8 



Q, > 15000 GeV 



10 

5 



1 1 1 1 1 1 1 


1 1 1 1 


id — 
















HI 


1 


1 1 1 1 1 1 1 1 1 



0.2 0.4 0.6 0.8 



M„> 180 GeV 




0.2 0.4 0.6 0.8 



Figure 3: Distributions of Me and ye for the selected NC DIS candidate events, (a) and (b) for 
2500 <Ql< 15000 GeV^ (c) and (d) for Ql > 15000 GeV^; distributions of (e) for ye > 0.4 and 
(f) ye for Me > 180 GeV; superimposed on the data points (• symbols) are histograms of the standard 
NC DIS expectation. 



The ye distribution of these high mass events is shown in Fig. 

Hence, no significant excess is seen in the mass spectrum for the 1997 data alone and the 
"clustering" around Me ~ 200 GeV is, overall, thus rendered less significant than that observed 
with 1994 to 1996 data only. 

We then consider the information. Figure ^ja shows the measured Ql distribution in 
comparison with the expectation from standard NC DIS processes. Also shown in Fig. ^ is 
the ratio of the observed Ql distribution to the NC DIS expectation. The errors resulting from 
the convolution of the systematic errors and the statistical error of the Monte Carlo sample are 
correlated for different Ql bins and are indicated in Fig. ^ as lines above and below unity 
joining the ±la errors evaluated at the centre of each bin. These errors are dominated by the 
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Figure 4: (a) Ql distribution of the selected NC DIS candidate events for the data (• symbols) and for 
standard NC DIS expectation (histogram); (b) ratio of the observed and expected (from NC DIS) number 
of events as a function of Ql; the lines above and below unity specify the ±1(7 levels determined using 
the combination of statistical and systematic errors of the DIS expectation. 



uncertainty in the electromagnetic energy scale of the calorimeter and vary between 7.7% at 
low Ql and 25% at the highest values of Q^. The NC DIS expectation agrees well with the data 
for Ql < 10000 GeV^ while at larger Ql, deviations are observed, with a slight deficit around 
15000 GeV^ and a number of observed events at > 15000 GeV^ in excess of the NC DIS 
expectation. For > 15000 GeV^, 22 events are observed while 14.1 ± 2.0 are expected from 
standard NC DIS. 



5.2 Charged Current Deep-Inelastic-like Signatures 
5.2.1 Event selection and kinematics 

The inclusive selection of CC DIS-like events requires: 

1. no 6=*= candidate with Et > 5 GeV found in the LAr calorimeter; 

2. the total missing transverse momentum PT,miss > 30 GeV. 

These cuts eliminate the photoproduction and NC DIS background. To deal with specific back- 
ground sources to CC DIS, it is required that there be no isolated track with > 10 GeV 
found within the angular range 10° < 9 < 145°. This reduces the remaining contamination 
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to < 0.3% from misidentified NC DIS events where the positron has been scattered through a 
crack of the calorimeter and also suppresses eventual background from single W boson, while 
causing negligible efficiency losses for the CC DIS selection. 

The Q^, y and M are calculated using the Jacquet-Blondel ansatz [ P3| ] by summing over all 
measured final state hadronic energy deposits using: 

This method will henceforth be called the hadron method (/i-method). 

In addition to the cuts (1) and (2) above, the analysis is restricted to the kinematic domain 
Q\ > 2500 GeV^ and < 0.9. The resolutions in both and Ql degrade with increasing y 
since both SMh/Mh and 5Q\/Q\ behave as 1/(1 — yh) for yh ~ 1. Hence the high yh domain 
is excluded. Throughout the remaining domain, the CC trigger efficiency is 96.5 ± 2%. For 
LQ u + q events, these selection criteria ensure typical efficiencies varying between ~ 32% 
and ~ 79% for LQ masses ranging in 75 to 250 GeV. 

Following this CC selection, 213 DIS event candidates are accepted in good agreement with 
the standard CC DIS expectation of 199.1 ± 11.5 events. 

5.2.2 Comparison with Standard Model expectation 

Figure |]a shows the two dimensional distribution of yh against for the CC candidates. 
Signal Monte Carlo events coming from the decay of a 200 GeV narrow scalar resonance into 
u + q are shown in the same plane in Fig. |]b, where the degradation of the resolution in 
at high yh is clearly visible. This resolution on the LQ mass is of about 10%. While the 
relative calibration of the hadronic scale is known at the 2% level (as discussed in section |^ 
and controlled comparing real and simulated NC DIS events), the energy scale procedure which 
relies on transverse momentum balance does not attempt to correct on average in DIS events 
the measured Mh to a "true" value. The measured Mh underestimates the resonance mass 
systematically by ~ 6% and this shift will be taken into account in deriving LQ results. 

Figures ^ and ^ show for the CC selection the measured Ql distribution in comparison 
with the standard CC DIS expectation. As can be seen in Fig. |]b, the systematic errors are rel- 
atively large and dominated by the uncertainty on the hadronic energy scale of the calorimeter. 
In the kinematic region Q\ > 15000 GeV^, there are Nobs = 7 observed events compared with 
an expectation of 4.84 ± 1.42 from standard CC DIS. 

5.3 High Pt Muon Signatures 

For LQs possessing a coupling to a second generation lepton, leading to fi + q final states, 
we start from a CC-like selection in the kinematic domain Ql > 1000 GeV^ and PT,miss > 
25 GeV, and then require: 
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Figure 5: Kinematics in the t/h - Mh plane of (a) the selected CC DIS candidate events from HI data 
(two isocurves at = 2500 and 15000 GeV^ are plotted as full lines); (b) a scalar leptoquark resonance 
of mass Mlq = 200 GeV decaying into v + q, for a coupling A = 0.05. 



1. an identified isolated muon with transverse momentum > 10 GeV and within polar 
angular range 10° < 6^ < 145°, where the momentum and angle are determined by the 
associated track. There should be no other charged track linked to the primary interac- 
tion vertex within the pseudorapidity-azimuthal isolation cone centred on the /i track of 
opening ^{/^ri^f + (A0^)2 = 0.5; 

2. at least one jet found in the angular range 7° < 9 < 145°, with a transverse momentum 
Pt > 15 GeV, using the cone algorithm mentioned in |5.1.1| . 

The muon identification combines inner tracking and calorimetric information. Within the 
fi isolation cone, at least 1 GeV must be visible in the calorimeters. Restricting to the LAr 
calorimeter, this energy should be smaller than one third of the fi track momentum. Less than 
5 GeV should be seen in the conical envelope between 0.5 < a/ (Ar]^)^ + (A^^)^ < 1.0. The 
centroid of the energy deposits in the calorimeters within the /i isolation cone should not be 
in the LAr electromagnetic section. With these identification criteria, muons are found with a 
typical efficiency of ^ 85% over most of the angular range. 
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Figure 6: (a) Qf^ distribution of the selected CC DIS candidate events for the data (• symbols) and for 
standard CC DIS expectation (histogram); (b) ratio of the observed and expected (from CC DIS) number 
of events as a function of Q^; the lines above and below unity specify the ±la levels determined using 
the combination of statistical and systematic errors of the DIS expectation. 



Only four fi + jet events are observed after applying these basic requirements. These four 
"outstanding" events are amongst^ the high Pt lepton events discussed in [p4|]. With these 
selection criteria 0.60 ± 0.10 events are expected from SM processes, coming mainly from W 
production and inelastic photon-photon interactions 77 /i+yU~. 

For the 2 — > 2 body LQ induced processes, where the final state consists of only the muon 
and the scattered quark, the energy-momentum conservation relates the polar angle 9i of the 
muon to the measurement of the hadronic final state, by : 

The selection of /i + jet candidates induced by LQ decay or w-channel exchange requires that : 

(i) the polar angle 9fj_ of the high Pt isolated track agrees with 9i within 30°; 

(ii) the track and the hadronic flow are back-to-back (A0^_/i > 170°) in the plane transverse 
to the beam axis. 



As seen in Fig. 0, these criteria induce a minute efficiency loss for a LQ signal simulation. 
The actual loss was estimated to be ^ 5% using real NC-like data by searching for the charged 



The event labelled /i3 in [Q], in whi ch bot h an isolated muon and the scattered positron are identified, does 
not fulfil the first requirement (see section 5.2.1 1 we apply to select CC-like events. 
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Figure 7: (a) Azimuthal opening between the track associated to the high Pt muon in e + p — > /i + jef 
events and the hadronic energy flow; (b) polar angle opening between the muon track and the final state 
lepton angle predicted from the hadronic energy flow for a 2 ^ 2 body process. Arrows indicate the A0 



and A6 values for the "outstanding" events reported in [44] and labelled as therein. Histograms represent 
the expected distributions of these variables for eq — > iJ,q' LQ Monte Carlo events. 



track angles associated with the positron candidate as deduced using the hadronic energy flow. 
The above requirements lead to typical efficiencies to select LQ ^ + q events which vary 
between ~ 30% and ~ 60% for LQ masses ranging in 75-250 GeV. For LQ masses far above 
the kinematic limit {Mlq ^ ^/s^, the distribution of the polar angle of the final state lepton 
depends on the LQ spin and fermion number (this latter affecting the relative contributions of 
the s- and w-channels). With the above selection criteria the efficiency to select fi + jet events 
induced by very heavy LQs varies between ~ 20% and ~ 33% depending on the LQ quantum 
numbers, but independently on the LQ mass. 

We observe no candidate satisfying the LQ fi + jet selection while 0. 12 ± 0.05 are expected 
from SM processes (mainly from inelastic photon-photon interactions). In particular, the e^p — > 
/i+X outstanding events discussed in [Q| fail significantly the kinematic constraints for LQ 
induced eq fxq' processes as can be seen in Fig. 0. 
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5.4 High Pt Tau Signatures 



For LQs possessing a coupling involving a third generation lepton leading to r + g final 
states, the analysis is restricted to hadronic decays of the r [|. The identification of "pencil-like" 
jets induced by hadronic decays of r requires that the jet invariant mass satisfies Mj^t < 7 GeV 
and has a low multiplicity, namely 1 < N tracks < 3, N tracks being the number of vertex fitted 
tracks with > 0.15 GeV in the jet identification cone. The jet mass Mj^t is here calculated 
as the invariant mass of all energy deposits associated to the jet, each of those being treated as 
a massless object. 

Inclusive r + jet signatures, where the r lepton decays into hadrons, are selected by requiring 
that : 

1. no candidate with Et > 5 GeV is found in the LAr calorimeter; 

2. two jets are found in the angular range 7° <9 < 145° using the cone algorithm mentioned 
in |5.1.1| with a transverse momentum Pt > 30 GeV; one of these jets must satisfy the 
"loose" r-jet identification criteria described above; 

3. there is at most a small amount of energy deposited in the backward calorimeter, Eback < 
7.5 GeV; 

4. the impact point of the "leading track" associated to the r-jet candidate at the inner surface 
of the LAr calorimeter must be at least 2° apart in azimuth from each of the eight </> 
cracks of the LAr. The "leading track" is, among the vertex fitted tracks found in the jet 
identification cone, the one which has the highest momentum projected on the jet axis; 

5. the leading track associated to the r-jet candidate carries a large enough fraction Etrack / Ejet > 
10% of the jet energy Ejet', 

6. the fraction of the r-jet candidate energy carried by the leading track and the energy 
deposition fraction fem in the LAr electromagnetic section of the r-jet are such that 

fem ~l~ Etrack / Ejet ^ 1'5 j 

7. the energy deposits of the jet should present an important longitudinal dispersion Rl = 
a/ {P) — {ly > 7 cm where, for each deposit, I is the distance from the impact point at 
the calorimeter surface along the jet axis; 

8. at most two vertex fitted tracks are found in the azimuthal hemisphere containing the r-jet 
candidate. 

Cuts (1) to (3) ensure a preselection of r -i- jet candidate events at high Pt; low NC DIS is 
suppressed by cut (3). At this stage, 375 events are selected in the data while 365.7 are expected 
from NC DIS and 7p processes. Figure ^ shows, for these selected candidates, the distribution 
of the polar angle Or-jet of the r-jet candidate. Photoproduction background, where a low 

^ The r+ jj^i'^pT channel is covered implicitly by the above n + q search since the t decay products are 
strongly boosted in the t direction; the r+ e^v^Vr channel is not covered here. 
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Figure 8: Distributions of (a) the polar angle Or- jet of the r-jet candidate; (b) the transverse momen- 
tum of the T-jet candidate; (c) the total missing transverse momentum; (d) the fraction of the energy of 
the T-jet candidate carried by the leading track and (e) the number of tracks in the azimuthal hemisphere 
containing the r candidate. Symbols correspond to data events and histograms to SM simulation. 



multiplicity jet can fake a r-jet, is seen to contribute mainly at low 9T—jet in contrast to NC DIS 
background arising when the scattered positron has not been identified. NC DIS contamination 
at low values of 9r-jet is due to events where the positron has been scattered at large angle 
through a crack of the calorimeter, and where the "current jet" has been identified as a high Pt 
r-jet. Figures and c show the distributions of the transverse momentum of the r-jet and of 
the whole final state. Figure ^ shows the ratio of the energy of the leading track associated with 
the r-jet candidate, to the r-jet energy. The observed distribution of this fraction Etrack/ Ejet is 
shifted by 9% compared to the one expected from the simulation. This shift is taken into account 
as a systematic error on the quantity Etrack/ Ejet, propagated when estimating the uncertainty 
on the SM expectation. On Fig. |8|e the number of tracks found in the azimuthal hemisphere 
containing the r-jet candidate is shown to be well described by the simulation. 

NC DIS and 7p backgrounds are further reduced by requirements (4) to (8). Cut (4) avoids 
regions close to (p cracks of the calorimeter and hence suppresses NC DIS background. Cut 
(5) efficiently reduces 7p background where a low multiplicity jet could fake a r-jet. Cut (6) 
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Figure 9: Correlation between the 
electromagnetic fraction of the r-jet 
candidate and the fraction of the jet en- 
ergy carried by the leading track for (a) 
NC DIS and (b) LQ ^ r+q simulated 
events; (c) shows the sum of these 
variables and the cut applied. Events 
where the track associated to the r-jet 
candidate is at less than 2° in azimuth 
from a (p crack of the calorimeter have 
been excluded. 
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ensures a powerful suppression of the remaining NC DIS background events where the scattered 
positron has not been identified. It exploits the fact that the electromagnetic fraction of a r-jet is 
high mainly when the r lepton decays to nn^ + X, which generally implies a small value for the 
ratio Etrack/ Ejet- On the contrary, for NC DIS events where the positron has been misidentified 
because its shower has a substantial leakage into the hadronic section of the calorimeter, this 
ratio Etrack I 'El jet is expected to peak at one. Figure ^ illustrates this property of r jets and shows 
how cut (6) discriminates the searched signal from the DIS background. Remaining NC DIS 
events are further removed by cut (7). The more stringent requirement on the track multiplicity 
imposed by cut (8) further reduces the remaining 7p contamination. This cut takes into account 
the fact that two reconstructed tracks can be associated to a single charged particle, especially 
when scattered in the forward region. 

The efficiency to identify r leptons decaying hadronically using the above criteria is ~ 25%. 

We observe 21 events satisfying the above inclusive r -i- jet requirements, which agrees well 
with the mean expectation of 25.4 ± 4.3 events coming from NC DIS and 7p processes. 

For the r + X channel, additional cuts relevant for the specific LQ search are applied : 

(i) a total transverse momentum PT-miss > 10 GeV ; 

(ii) the r-jet candidate is at Acpr-h > 160° from the total hadronic flow. 
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Both these cuts exploit the fact that neutrino(s) emerging from the r decay are collimated with 
the r-jet direction, and that the other jet in LQ induced r+jet events is formed by the fragmenta- 
tion of a quark and thus should not contribute to the missing transverse momentum. Moreover, 
cut (i) efficiently reduces the remaining NC DIS contamination. The correlation between the 
Premiss and the /^(pr-h for the 21 data events which satisfy the criteria (1) to (8) listed above is 
shown in Fig. |m together with the expected distributions of /^(pr-h and PT,miss for LQ induced 
r+jet Monte Carlo events. 
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Figure 10: Distributions of (a) 
the azimuthal angular opening ^(p^-^^ 
between the r-jet candidate and the 
hadronic energy flow; (b) the missing 
transverse momentum PT,miss- White 
and hatched histograms correspond to 
the SM expectation and the LQ in- 
duced r-f-jet signal respectively, (c) 
Correlation of these two variables for 
the data events satisfying cuts (1) to 
(8). The full lines show the cuts ap- 
plied, which retain the upper-right part 
of this plane. 



The final selection efficiency on the LQ t + q signal is ~ 10% for LQ masses of 100 GeV 
and reaches a plateau at ~ 25% above 200 GeV. The small efficiency at low masses is mainly 
due to the requirement of two high Pt jets. For LQ masses far above the kinematic limit, this 
efficiency varies between ~ 8% and ^12% depending on the LQ spin and fermion number. 

We observe no candidate satisfying the additional cuts designed specifically for LQ in- 
duced processes involving third generation leptons for an expectation of 0.77 ± 0.30 event from 
misidentified electrons in NC DIS processes. 

It should be noted that the absence of r -i- jet candidates accompanied by a large PT,miss, 
despite the relatively loose r-jet requirements, makes it rather unlikely for the yu's in the e'^p 
/i+X events [ p4| ] to originate from a jet fluctuating in one single particle. Conversely, these 
muon events fail significantly the kinematic requirements of a two body decay LQ ^ q + r 
followed by a r /i+z/^z/^ decay since the final state /i+ would appear boosted in the r direction 
and the A9-A(j) restrictions of Fig. would still apply. 
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6 Constraints on First Generation Leptoquarks 



6.1 The Leptoquark Specific Angular Cut 

We first consider LQs possessing couplings to e — g or e — g pairs only. 

In order to enhance the significance of a possible LQ signal over the NC DIS background 
remaining after the selection requirements described in section the specific angular distri- 
butions for LQ induced processes are exploited. For scalar or vector LQs with either F = or 
I F 1= 2, distinct mass dependent lower i/g > Ucut cuts have been optimized using Monte Carlo 
generator programs to maximize the signal significance. This has been achieved by finding the 
best compromise between the efficiency loss on the LQ signal and the important background 
reduction. Thus, with increasing LQ mass, the Ucut decreases together with the NC DIS expecta- 
tion. However, for LQ masses close to the kinematic limit, the mass spectrum is highly distorted 
towards low values (as explained in section ^. Thus, a higher Ucut will be needed to enhance 
the signal significance of such a high mass LQ observed at -C M^q where a larger NC DIS 
background is expected. The evolution of Ucut as a function of the LQ mass is represented in 
Fig. PT]a, for the example of F = and | F |= 2 scalar LQs. For a scalar F = LQ, the Ucut 
monotonously decreases from ~ 0.6 around 60 GeV to ~ 0.2 around 260 GeV, and then rises 
up to ~ 0.5 at 300 GeV. The behaviour is similar for a scalar | F |= 2 LQ, but the effect of 
the distortion of the mass spectrum is significant already at ~ 200 GeV. A similar description 
holds for vector LQs but smaller values of the ycut are obtained, because of the [1 — yY shape 
of their y spectra. 



6.2 The Mass Spectra 

The comparison of the measured mass spectrum with SM predictions is shown in Fig. |T2|a for 
the NC DIS analysis. The measured and expected NC spectra are seen before and after applying 
the mass dependent ye cut relevant for a F = scalar LQ. After this y^ cut, we observe 310 
events in the mass range Me > 62.5 GeV while 301.2 ± 22.7 are expected. The observed mass 
spectrum is seen to be well described by the SM expectation, with nevertheless a slight excess 
in the mass range 200 GeV ± AM/2 with AM = 25 GeV, due to the same 8 events already 
discussed in section |5.1.2| (the cut optimal for F = LQs is y^ut — 0.4 for M ~ 200 GeV). 
For CC DIS-like final states, the mass distribution of the 213 events fulfilling the requirements 



described in section 5.2.1 is seen in Fig. 12b to be in good agreement with the SM prediction. 



6.3 The Limits Derivation 

Assuming that the observed excess of events in the NC DIS-like channel is due to a statistical 
fluctuation (or either, formally, that the event sample contains at most two components, an 
unknown LQ signal and a known expectation from NC DIS), models containing first generation 
leptoquarks can be constrained. 

An upper limit Num on the number of events coming from leptoquark induced processes 
can be obtained assuming Poisson distributions for the SM background expectation and for the 
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Figure 1 1 : Lower cut applied (a) on ye and (b) on Me to enhance the signal significance for a F = 
(full line) and a | F | = 2 (dashed line) scalar leptoquark, for NC DIS-like final states; contribution aiQ 
of the LQ induced e+p —>■ e + q + X processes (s- and u-channels summed) and of their interference 
Gint with SM DIS for (c) Sq^r (| F |= 2) and (d) S'i/s^l {F = 0) leptoquark, with A = 0.5. In (c), 
the contribution of the u-channel LQ exchange alone is also represented (dotted line). Each contribution 
has been integrated over the phase space allowed by the mass dependent Mg and cuts applied in the 
analysis. 

signal. For each contributing channel, we use the numbers of observed and expected events 
within a mass bin [Mmin-, Mmax] of variable width, adapted to the expected mass resolution and 
measured mass values for a given true LQ mass, and which slides over the accessible mass 
range. For example, only NC DIS (respectively CC DIS) candidates G [186; 204] GeV 
(Mh e [170; 214] GeV) will be used to constrain a 200 GeV LQ undergoing a NC (CC) DIS- 
like decay. For high LQ masses, the mass bin becomes very large because of the distortion of 
the mass spectrum mentioned in section |[ The evolution of Mmin as a function of the LQ mass 
is represented in Fig. |TT]b for F = and | F |= 2 scalar LQs in the NC-like channel. For LQs 
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Figure 12: Mass spectra for (a) NC DIS-like and (b) CC DIS-like final states for data (symbols) and 
DIS expectation (histograms). In (a), the NC DIS-like comparison is shown before (open triangles, white 
histogram) and after (closed dots, hatched histogram) a y cut designed to maximize the significance of 
an eventual scalar LQ signal. The greyed boxes indicate the ±lcr band combining the statistical and 
systematic errors of the NC and CC DIS expectations. 

undergoing NC DIS-like decays, typical signal detection efficiency including the optimized y 
and mass cuts is found to vary between 24% at 75 GeV, ~ 35% around 200 GeV and 48% 
at 250 GeV. In the CC DIS-like decay channel, this efficiency ranges between ~ 19% for a 
75 GeV leptoquark, ~ 39% at 100 GeV and reaches 53% at 200 GeV. 

In section |64| , we will consider the BRW model where the only free parameter is the Yukawa 
coupling A. Only NC DIS-like data will be used to derive Num, which can then be translated 
into an upper limit on aiq + aint, ctlq being the part of the e+p ^ e + q + X cross-section 
induced by the LQ s- and w-channel processes, and a^nt the interference term with DIS, both 
integrated over the phase space allowed by the mass dependent cuts applied in the analysis on 
Me and ye. The evolution of these integrated contributions aiQ and a^t as a function of the 
LQ mass is shown in Fig. |lT]c for the Sq^r (| F |= 2) scalar leptoquark, as the dashed and 
dash-dotted lines respectively, and in Fig. [TTI d for the 5'i/2,l (F = 0), with emphasis on the 
high mass domain where the interference contribution becomes visible. In the former case the 
interference between LQ induced processes and SM boson exchange is constructive, while it is 
destructive for the latter. The cross-sections have here been calculated for a fixed value A = 0.5, 
which is typical for the experimental sensitivity in the displayed mass range. In Fig. [TTI c the 
contribution of the w-channel Sq^r exchange alone is also shown as the dotted line (for the 
case depicted in Fig. [TT|d, cr^ is below 10^'* pb and is not represented in the figure). As was 
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mentioned in section 3, the mass dependent cuts applied on Me and Ue considerably reduce the 
contributions of the interference and of virtual exchange (e.g. by a factor (9(10) for a Sq^r LQ 
at M iQ = 250 GeV). As can be seen in Fig. pT] c and d, the interference is negligible for F = 
leptoquarks and masses up to ~ 275 GeV, but plays an important role for \ F \= 2 LQs as soon 
as Mlq ~ 220 GeV. Moreover the n-channel contribution is always negligible. 

In the mass domain where the interference between standard DIS and LQ induced processes 
can be neglected, we will move away from the BRW model and consider a more general case 
where the branching = P{LQ eq) is not determined by A only. Taking (3^ and jSy = 
P{LQ uq) as free parameters, NC and CC DIS-like data can be combined to derive Num, 
which can then be translated into an upper limit on the signal cross-section o-^q = ag + (Tu 



and thus on the Yukawa coupling A. This will be done in section p3| . The signal cross-section 
in the mass domain considered being largely dominated by the s-channel resonant production 
as mentioned above, an upper limit on ag x /3e can also be derived. Fixing A, mass dependent 
upper limits on the branching ratio Pe can then be obtained, as will be done in section |6^ using 
NC DIS data only. In these two cases, we make use of the signal detection efficiencies given 
above to translate into an upper limit on the signal cross-section. 

The procedure which folds in, channel per channel, the statistical and systematic errors is 



described in detail in [PTI] 



6.4 Mass Dependent Limits on the Yukawa Coupling in the BRW Model 

We first establish constraints on the BRW model described in section ^ taking into account all 
LQ induced contributions, but restricting the analysis to NC DIS-like processes. For the decay 
of resonantly produced LQs, the values of P^, are specified and given in table |^. As discussed 
in section ^, the constraints can be extended beyond the kinematic limit by profiting from the 
tail expected in the s-channel towards low masses, and by properly taking into account the 
interference with SM boson induced processes. 

In the very high mass domain, the interference 0-^4 of LQ induced processes with NC DIS 
generally dominates the LQ cross-section aig- Instead of a^q alone, we are thus directly 
sensitive to (Jsum = cr^g + crint, which, for those LQ species which interfere destructively with 
DIS, could be negative when integrated over the whole phase space, but remains positive within 
the kinematic cuts applied. We thus proceed the following way : 

• For a given LQ mass Mlq, the numbers of observed and expected events within the 
optimized cuts M^in < < M^ax and j/e > Ucut are used to set an upper limit 
on the number of signal events. For LQ masses above the kinematic limit, the optimized 
Me and ye cuts are nearly independent of M^q and close to those displayed in Fig. [TT|a,b 
for Mlq = 300 GeV. A first estimate Aq of the upper limit on A is then obtained by 
solving : Num = sum,cuts{>^) ■ crsum,cuts is calculated analytically by integrating over 
the phase space allowed by the cuts the squared amplitude for the LQ process and its 
interference with NC DIS. The cuts applied, by reducing the NC DIS contribution, ensure 
that cisum,cuts is positive for reasonable A values. 
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The event generator LEGO is then used to produce LQ events at (M^g, Aq), from which 
we get the acceptance correction factor A defined as the ratio of the number of events 
satisfying the NC DIS selection cuts as well as the Me and cuts, to the number of 
events generated within these cuts. 

To take into account next-to-leading order QCD corrections on the LQ production cross- 
section, we calculate the convolution K* of the K-factor given in for the NWA, 
with the LQ Breit-Wigner distribution corresponding to {Mlq, Aq). The resulting func- 
tion K*{Mlq) will be used henceforth to calculate the signal cross-sections. The NLO 
corrections lead to a sizeable enhancement of the cross-section as will be seen below. 



The final rejected coupling is then obtained by solving the equation Num 

K* X 0'sum,cuts{^) ■ 



£x Ax 



The resulting rejection limits at 95% confidence level (CL) are shown in Fig. [T3| up to M^q = 
400 GeV for | F |= or 2 scalars or vectors. Constraints for even larger masses where one 
is probing "strong" (i.e. A > 1) coupling values in a contact interaction will be discussed in a 
separate paper. The case of | F | = LQs, which can be produced via fusion between the e"*" 
and an incident valence quark, naturally offers the best sensitivity. 




10 



10 



100 150 200 250 300 350 400 

Mlq (GeV) 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
F = 


1 1 1 1 1 1 1 1 1 




(b) : 




fl'2, L 




■■■■ |l/2,L 
*1/2, R 


■v-r-T-;-;" , , 1 , , , , 1 , , , , 1 , , 


HI 



100 150 200 250 300 350 400 

Mlq (GeV) 



10 



10 



10 



10 



-2 



^ |F|=2 






(c) : 




— V 

* 1 / 2, L 

^ 1 /2, L 
... w 

Hln 



100 150 200 250 300 350 400 

Mlq (GeV) 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 
F = 


1 1 1 1 1 1 1 1 1 
MMJMJMi 






(d) ; 


-\j}mimM>C---'' .... />■ 

'-:":'\ . . . . \ . . . 


V " " Yo, R 
- V, L 


HI 



100 150 200 250 300 350 400 

Mlq (GeV) 



Figure 13: Exclusion limits at 95% CL on the Yukawa coupling A as a function of the leptoquark mass 
for (a) I F I = 2 and (b) F = scalar, (c) | F | = 2 and (d) F = vector leptoquarks described by the 
BRW model. Domains above the curves are excluded. 



For LQs having F = 0, these limits represent an improvement by a factor ~ 3 compared 
to HI published results [pT|]. For a Yukawa coupling of the electromagnetic strength A^/47r = 
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aEM (i-e. A ~ 0.3), such scalar (vector) LQs are excluded at 95% CL up to 275 GeV (284 GeV) 



by this analysis. As can be seen in Fig. [Tjb the Si/2,r is the scalar for which HERA has the 
highest sensitivity since both charge states can be produced via a fusion e+w or e^d. On the 
contrary, only an e^u (e+rf) fusion is possible for the production of 51/2, l (5*1/2, l), for which 
the cross-sections are thus smaller. Note also that, due to the more favorable parton density, 
a higher cross-section (and thus a better rejection limit) is expected for 5i/2,l than for 5'i/2,l- 
However at very high masses, the 5i/2,l interferes destructively with NC DIS, on the contrary 
to the 5i/2,L, so that the limits become similar. 

For I F 1= 2 leptoquarks, the improvement compared to our previous published results is 
less substantial since e^p data collected in 1994 had been also taken into account in [[TTl]. These 
LQs will be best probed with e'p data which are being collected since 1998. 

For coupling values equal to the obtained limits, NLO QCD corrections enhance the produc- 
tion cross-section of a F = LQ by ~ 10% for Mlq = 100 GeV to ~ 30% at 250 GeV. For 
higher masses this enhancement decreases because of the distortion of the LQ mass spectrum. 
For I F 1= 2 LQ, this enhancement remains below ~ 20%. 

5i,i and 5o,l leptoquarks being allowed to undergo CC DIS-like decays with a branching 
ratio Pi, = 0.5, the combination of NC and CC DIS analysis is expected to better constrain these 
leptoquarks. The result of this combination is given for the Sq ^ by the curve Pe = Pu = 0.5 in 
Fig. It can be seen that combining the two contributing channels improves the sensitivity 
on 5o,L up to that achieved on Sq^r, recalled on Fig. [T^a by the curve /5e = 1. 

The DELPHI [ff5]] experiment at LEP recently performed a direct search for single lep- 



toquarks using data accumulated at e+e centre of mass energies y/Se+e- up to 183 GeV. 
Constraints relevant for some LQ species have also been obtained indirectly by ALEPH [HQ], 



OPAL pT]] and L3 Collaborations from measurements of hadronic cross-sections and 



asymmetries at ^/s^+^ = 130 to 183 GeV. For a Yukawa coupling of the electromagnetic 
strength, the direct search via single production results in a limit of 171 GeV which is not yet 
competitive with this HI analysis. For similar coupling values, the indirect search at LEP is 
yet only competitive for the 5o,l, 5i,l, Vq^l, ^o,/? and Vi^l, for which limits up to 240 GeV and 
470 GeV have been obtained in the scalar and vector case respectively. However the sensitivity 
of the indirect search drops quickly with the Yukawa coupling, and for A = 0.1 bounds lie 
below 100 GeV for all LQ species 



In constrast to an ep collider, the limits derived at the TeVatron where LQs are mostly 
produced in pair via the strong interaction, with larger cross-sections expected for vector than 
for scalar leptoquarks, are essentially independent of the Yukawa coupling. Recent results have 
been published for scalar leptoquarks by D$ [^] and CDF [ |5m Collaborations. In the BRW 
model, a combination of these results [ pT] ] excludes scalar leptoquark masses up to 242 GeV. 
Comparisons of our results with limits obtained at the TeVatron collider in more general models 



will be presented in sections and ^ 



6.5 Mass Dependent Limits on the Yukawa Coupling in Generic Models 

Moving away from the BRW model, we now consider leptoquarks for which the branching ra- 
tios Pe and Pi, in NC and CC DIS-like decays are free parameters. As an example, the 95% CL 
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exclusion contour for a scalar LQ with \ F \= 2 decaying with (5i, = 90% in uq and Pe = 10% 



in eq is shown in Fig. [IJa as the dash-dotted line. In this case, limits are shown only in the mass 
domain where the interference of LQ processes with DIS can be neglected. The gain obtained 
by the combination of both channels can be seen when comparing this contour with the greyed 
domain showing the sensitivity achieved using only NC DIS-like channel. For Pi, = 90% and 
Pe = 10%, scalar LQ masses below 200 GeV are excluded at 95% CL by our analysis, for 
Yukawa couplings of the electromagnetic strength. This extends far beyond the domain ex- 



cluded by TeVatron experiments [^, which for such small values of P^, only exclude scalar 
leptoquark masses below ~ 110 GeV. At HERA the sensitivity will be considerably enhanced 
with an electron incident beam where these LQs can be produced via a fusion with a valence 
u quark. Hence a large discovery potential is opened for HERA, for high mass leptoquarks 
decaying in uq with a high branching ratio. 
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Figure 14: Exclusion limits at 95% CL on the Yukawa coupling A as a function of the leptoquark mass 
for (a) scalar leptoquarks produced via e^u fusion and (b) vector leptoquarks produced via e^d fusion. 
Different hypotheses for the branching ratios into eq, vq are considered. Domains above the curves are 
excluded. 



Similar results are given for a vector LQ coupling to e^d in Fig. [T^. For the above values of 
(Pe, Pu) and of the Yukawa coupling, the excluded mass domain extends in this case to 250 GeV. 



6.6 Mass Dependent Limits on f3{LQ eq) 



We consider here leptoquarks which undergo NC DIS-like decays with a branching ratio P^, and 
do not make any assumption on the other possible decay modes of the LQ. For a fixed value of 
the Yukawa coupling A upper limits on LQ production cross-section are translated in terms of 
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mass dependent limits on the branching Exclusion limits at 95% CL are shown in Fig. [T3|a,b 



as greyed areas for scalar LQs produced via a e^d or e+w fusion respectively. On Fig. [T^a, the 



CQ. 



Scalar LQ coupling to + d 




1 



+ 

0) 

T 
o 



75 100 125 150 175 200 225 250 275 

Mlq (GeV) 

Scalar LQ coupling to + u 



0.8 
0.6 
0.4 
0.2 

0. 



DO combined Limit 



□ HI limit ( X=0.1 ) 

□ HI limit ( X=0.05 ) 



(NC data only) 




100 125 150 175 200 225 250 275 

Mlq (GeV) 



Figure 15: Mass dependent exclusion limits at 95% CL on the branching ratio f3{LQ — > eq) for scalar 
leptoquarks produced by (a) e~^d and (b) e~^u fusion. Two exclusion regions corresponding to A = 0.1 
and A = 0.05 are represented as greyed areas. For A = 0.1, the error bands on the exclusion limits in (a) 
and (b) illustrate the sensitivity to d and u quark densities respectively (see text). The D$ limit is also 
shown as hatched region. 



central limit contour for A = 0.1 has been obtained assuming an uncertainty on the d quark 
density distribution which varies with x as mentioned in section ^. Curves above and below this 
contour define the error band obtained when the lack of knowledge on the proton structure is 
directly applied on the parton densities, instead of entering as a systematic uncertainty : the d 
quark density is enhanced/lowered by a factor varying linearly between 7% at low masses and 
up to 30% at 250 GeV. For LQs coupling to e+w (Fig. p3|b), a constant scaling factor of ±7% 
has been applied on the quite well-known u density, resulting in a much narrower error band. 

Despite the small A values considered, domains excluded by this analysis extend beyond 
the region covered by the experiment at the TeVatron also in the less favourable case 
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of LQs coupling to e~^d. This is especially the case for small values of /?e- For example, for 
Pe = 10%, this analysis rules out masses below ~ 240 GeV (~ 200 GeV) if A = 0.1 (A = 0.05) 
independently of other possible decay modes of the leptoquark, as can be seen in Fig. p3|b. This 
limit extends up to 260 GeV for an electromagnetic strength of the Yukawa coupling A as will 
be shown by the greyed domain in Fig. [T^. Hence, a yet unexplored region in the high mass - 
low Pe domain is probed by this analysis. 



7 Constraints on Couplings to Mixed Generations 

In this section, we consider LFV LQs which couple to both the electron and a second or third 
generation lepton. 

The case of low mass (M < y/s^) LFV LQs is first addressed in section [7?!] . The analysis is 
there restricted to LQs possessing a coupling An, allowing the LQ to be produced or exchanged 
between the incident lepton and a valence quark coming from the proton. Moreover we do 
not consider e ^ fi transitions induced by a low mass LQ, such processes being strongly 
constrained by low energy experiments as will be seen in section where high mass (M > 
y/s^) LFV LQs are considered. In this latter case, the study will be extended to any Ai,;, and 
both and X^j. 



7.1 Low Mass (M < ^/s^) LFV Leptoquarks 

We now consider LFV LQs possessing a coupling An to first generation lepton-quark pairs as 
well as a coupling X^j with leptons of the third generation. 

Mass dependent exclusion limits on An at 95% CL are shown in Fig. |T6|a when fixing 
the relative branching fractions Pe and Pr into e + jet and r + jet final states. A generic 
scalar leptoquark coupling to e"*" + m pairs (such as the 5'i/2,l in the BRW model) has been 
considered for three different sets of (P^,, Pr)- Here both e+ + jet and r+ + jet channels are 
combined. This latter channel is background free but the former benefits from a higher selection 
efficiency, such that finally both provide a comparable sensitivity to the signal. Thus, as soon 
as Pe + Pt approaches 1, the limits derived are very similar to those obtained for Pe = I. 
Assuming Pe = 10% and Pr = 90%, masses below 275 GeV are excluded at 95% CL for an 
electromagnetic strength of the An coupling. Such limits are especially interesting since, for 
small Pe and high P^-, the mass domain excluded by the TeVatron does not extend very far. 
The CDF experiment has performed a search for third generation LQs looking at rrbb final 



states [p2|], and excluded scalar LQs with masses below 99 GeV if P{LQ rb) = 1. A 
complementary search has been carried out by [|53[], where the analysis of uvbb final states 
leads to a lower mass limit of 94 GeV for P{LQ ub) = 1. HERA thus appears to provide 
access to an unexplored domain for LQs decaying with a small branching ratio in eg and a high 
branching ratio in rq. 

An alternative representation of our results is given in Fig. in the plane X-sj against the 
LQ mass, for different fixed values of An. We consider here a scalar LQ formed via e'^u fusion 
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Figure 16: (a) Mass dependent exclusion limits at 95% CL on the Yukawa coupling An, for scalar 
leptoquarks produced by e+n fusion. Different hypotheses for the branching ratios into eq, rq are con- 
sidered. Domains above the curves are excluded, (b) Exclusion domains in the plane (j = 1,2) 
against the leptoquark mass for several fixed values of An (greyed areas). A mass dependent indirect 
limit on A31 is also represented by the dashed line. 



(carrying the quantum numbers of the 5'i/2,l in the BRW model) such that only couplings Asj 
with j = 1, 2 are relevant and make the additional assumption that (3^ + = 1. The r + jet 
final states analysis provides a sensitivity on A31 so long as the LQ is light enough to have a 
substantial production cross-section via An, as can be seen on the 95% CL exclusion domains 
shown in Fig. |T^. For both couplings An and X^j of the electromagnetic strength, such LQs 
lighter than 270 GeV are excluded at 95% CL. This limit still reaches 237 GeV for LQs formed 
by e^d fusion (not shown). A similar lepton flavour violation analysis has been carried out by 
ZEUS Collaboration [ |T3| ] using an integrated luminosity of ~ 3 pb For An = As-, = 0.03, 
the analysis presented here extends their excluded mass range by ~ 65 GeV. Also shown in 
Fig. is the best indirect limit [ |22| ] on A31 in the case An = 0.03. This indirect constraint 



comes from the upper limit on the branching ratio P{t ^ vr^e) which could be affected by the 
process r+ — > m + LQ followed by the e'^ + u decay of the virtual leptoquark. Here the most 



recent upper limit [ ]54l ] on /3(r — > vr^e) has been used to update the bound derived in [p2|]. The 
HI direct limit improves this indirect constraint by typically one order of magnitude. No low 
energy process constrains the coupling products involving A32. In this case, HI covers a yet 
unexplored domain. 

It should be noted however that for most other LFV LQ species (except namely Vq^r, Sq^r 
and Vi/2,l) the coupling products involving A32 and A33 can be constrained by low energy ex- 
periments, in particular by r ^ K^e, B reX, Vub measurements ox K ^ -nvv W^- This 



latter process yields the most severe bound, relevant for An x A32, but which for F = LQs 
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only applies to those possessing the quantum numbers of the Vi^l in the BRW model. It has 
been checked (not shown here) that for all other F = LQs, the HI direct limits on these cou- 
pling products extend beyond the domain covered by low energy phenomena. More detailed 
comparisons of HERA sensitivity with indirect bounds will be discussed in the next section. 



7.2 High Mass (M > LFV Leptoquarks 

In this section, we make use of the fact that no /i + jet or r + jet candidate was found (with 
kinematic properties compatible with a 2 ^ 2 body process) to set constraints on very high mass 
LFV leptoquarks. For LQ masses well above the kinematic limit, the cross-section cr(e+ + p ^ 
C + + ^) depends only on {XuX^j / MIqY , with i and j indexing the generation of the 
quark coupling to LQ — e and LQ — /„ respectively, and where In = 1^ for n = 2 and In = t for 
n = 3. 

The 95% CL rejection limits are given for a e ^ /i transition in Fig. [T^ for F = and 
Fig. ITE| for I F 1= 2 leptoquarks. Results are obtained for all possible quarks involved, qi 



being the quark coupling to LQ — e and qj the one coupling to LQ — ji. The limits are given 
in units of 10"^ GeV^l For processes involving a b quark in the initial state (i = 3 or j = 
3), it has been checked that the correction to the cross-section due to the finite mass of the b 
remains below ~ 5% [^3p. Early HERA results were presented in a similar representation by the 
ZEUS Collaboration in [|13]]. The nomenclature defined in table ^ has been kept to distinguish 



all possible LQ quantum numbers. Also given for each entry in the depicted tables are the 



constraints from the indirect process which currently provides the most stringent bound p^]. 



Figures |T9| and ^ show similar limits for the e r transition. 



In Figs. ^7| - |20| , the bounds derived in [ P2| ] have been updated to take into account the latest 
results on suppressed or forbidden decays [|5^, and on the /i — e nuclear conversion [|5^]. 
Significant such updates concern in particular D and K decays into yue, B ^ fie as well as 
K Txvv, T Tie, r ^ 67 and /i ^ 67 [|5B|]. Previous HERA results [|T3]], now superseeded, 
are given in cases where there exists no known indirect constraint. 

Provided that the quarks involved do not both belong to the first generation, it is seen that in 
some cases HI limits supersede or come close to existing indirect bounds. For LQs coupling to 
muons, this concerns in particular the leptoquarks which can contribute to n ^ ej and D ^ fie 
processes. For LQs coupling to taus, this concerns more cases and in particular the leptoquarks 
which can contribute to r — > Ke, r — 67 and various forbidden B decays. The sensitivity 
improves over indirect constraints by up to an order of magnitude and some cases are covered 
uniquely by HERA experiments. 



8 Conclusions 

First generation leptoquarks (LQs) as well as leptoquarks possessing lepton flavor violating 
couplings have been searched at the HERA collider using 1994 to 1997 HI data in a mass range 
extending from 75 GeV to beyond the ep kinematic limit of ~ 300 GeV. 
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Figure 17: Rejection limits at 95% on XiiX2j/MlQ in units of 10~^ GeV~^ for F = leptoquarks, 
compared to constraints from indirect processes. The first column indicates the generations of the quarks 
qi and qj coupling respectively to LQ — e and LQ — ji. In each box, the process which provides a 
most [ p2| ] stringent indirect constraint is listed (first line) together with its exclusion limit (second line) 
and compared to the actual HI result (third line). Shadowed boxes emphasize where HERA limit is 
comparable to (within a factor of 2) or better than the indirect constraints. The open boxes marked with 
a * are cases which would involve a top quark. 



The search for leptoquarks which only couple to first generation fermions involved an anal- 
ysis of very high neutral (NC) and charged current (CC) deep-inelastic scattering data. The 
comparison of these data with Standard Model expectations has shown deviations in the 
spectrum at > iqooO GeV^ which are less significant than those previously observed in 1994 
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Figure 18: Rejection limits on AiiA2j /M|q in units of 10~^ GeV~^ for | F |= 2 leptoquarks, 
compared to constraints from indirect processes. The first column indicates the generations of the quarks 
qi and qj coupling respectively to LQ — e and LQ — ji. In each box, the process which provides a 
most [ p2| ] stringent indirect constraint is listed (first line) together with its exclusion limit (second line) 
and compared to the actual HI result (third line). Shadowed boxes emphasize where HERA limit is 
comparable to (within a factor of 2) or better than the indirect constraints. The superscripts Z indicate 
where earlier HERA [ pj| ] results already improved indirect bounds. The open boxes marked with a * are 
cases which would involve a top quark. 



to 1996. No significant clustering of events in excess of SM expectations has been found in the 
mass spectra for NC or CC-like events in the 1997 dataset alone. Exclusion domains for LQ 
masses and couplings have been derived. 
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Figure 19: Rejection limits on XuX^j / MIq in units of 10^^ GeV^^ for F = leptoquarks, compared 
to constraints from indirect processes. The first column indicates the generations of the quarks qi and qj 
coupling respectively to LQ—e and LQ—t. In each box, the process which provides a most [|^] stringent 
indirect constraint is listed (first line) together with its exclusion limit (second line) and compared to the 
actual HI result (third line). For the <S'i/2,L' which does not couple to the neutrino, limits on An x A32 
and on A12 x A31 derived in [ p2| ] from K tti/u have been replaced by the bounds obtained from 
r —>■ Ke. Shadowed boxes emphasize where HERA limit is comparable to (within a factor of 2) or 
better than the indirect constiaints. The superscripts Z indicate where earlier HERA [ 13 1 results already 
improved indirect bounds. The open boxes marked with a * are cases which would involve a top quark. 



For first generation leptoquarks of the Buchmiiller-Riickl-Wyler (BRW) effective model, 
masses up to 275 GeV (284 GeV) are excluded for scalars (vectors) with a Yukawa coupling 
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Figure 20: Rejection limits on XuX^j /MIq in units of 10"'^ GeV^^ for | F |= 2 leptoquarks, 
compared to constraints from indirect processes. The first column indicates the generations of the quarks 
Qi and Qj coupling respectively to LQ — e and LQ — r. In each box, the process which provides a 
most [ 22 1 stringent indirect constraint is listed (first line) together with its exclusion limit (second line) 
and compared to the actual HI result (third line). Shadowed boxes emphasize where HERA limit is 
comparable to (within a factor of 2) or better than the indirect constraints. The open boxes marked with 
a * are cases which would involve a top quark. 



of electromagnetic strength, A = y/AnoEM = 0.3. Constraints on the LQ couplings have been 
established for A < 1.0 for all LQ types for masses up to 400 GeV. For A = 0.3 but in generic 
models with arbitrarily small decay branching ratio /?e into NC-like final states, the exclusion 
domain extends to 260 GeV for Pe as small as 10%, far beyond the present reach of other 
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existing colliders. 



No event candidate has been found with either /x + jet or r + jet final states compatible 
with the production of a LQ with couplings mixing the first and second or third generation in 
the leptonic sector. The constraints derived on the Yukawa couplings extend for some LQ types 
and coupling products beyond the reach of other colliders as well as of low energy experiments. 
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